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Abstract 

The coupled motion is investigated for a mechanical system consisting of 
water and a body freely floating in it. Water occupies either a half-space or 
a layer of constant depth into which an infinitely long surface-piercing cylin¬ 
der is immersed, thus allowing us to study two-dimensional modes. Under the 
assumption that the motion is of small amplitude near equilibrium, a linear set¬ 
ting is applicable and for the time-harmonic oscillations it reduces to a spectral 
problem with the frequency of oscillations as the spectral parameter. It is essen¬ 
tial that one of the problem’s relations is linear with respect to the parameter, 
whereas two others are quadratic with respect to it. 

Within this framework, it is shown that the total energy of the water motion 
is finite and the equipartition of energy holds for the whole system. On this basis, 
it is proved that no wave modes can be trapped provided their frequencies exceed 
a bound depending on cylinder’s properties, whereas its geometry is subject to 
some restrictions and, in some cases, certain restrictions are imposed on the 
type of mode. 


1 Introduction 

This paper continues the rigorous study (initiated in [?]) of the coupled time-harmonic 
motion of the mechanical system which consists of water and a rigid body freely 
floating in it. The former is bounded from above by a free surface, whereas the 
latter is assumed to be an infinitely long cylinder which allows us to investigate two- 
dimensional modes orthogonal to its generators. The body is surface-piercing and no 
external forces acts on it (for example, due to constraints on its motion). The water 
domain is either infinitely deep or has a constant finite depth; the surface tension is 
neglected on the free surface of water whose motion is irrotational. The motion of 
the whole system is supposed to be small-amplitude near equilibrium which allows us 
to use a linear model. 

In the framework of the linear theory of water waves, the time-dependent problem 
describing the coupled motion of water and a freely floating surface-piercing rigid 
body was developed by pQ. However, his formulation was rather cumbersome, and so 
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during the second half of the 20th century the main efforts were devoted to various 
problems involving fixed bodies instead of freely floating ones (see the summarising 
monograph by El)- The cornerstone was laid by [51 himself who proved the first result 
guaranteeing the absence of trapped modes at all frequencies provided an immersed 
obstacle has a fixed position and is subject to a geometric restriction now usually 
referred to as John’s condition. In the two-dimensional case, it includes the following 
two requirements: (i) there is only one surface-piercing cylinder in the set of cylinders 
forming the obstacle; (ii) the whole obstacle is confined within the strip between two 
vertical lines through the points, where the surface-piercing contour intersects the free 
surface of water, the part of bottom (when the depth is finite) is horizontal outside 
of this strip. 

EH demonstrated that if condition (i) holds, then condition (ii) can be replaced 
by a weaker one. Namely, if the depth is infinite, then the whole obstacle must be 
confined to the angular domain between the lines inclined at 7r/4 to the vertical and 
going through the two points, where the surface-piercing contour intersects the free 
surface. If the depth is finite, then it is required that the whole obstacle is confined 
to a smaller angular domain between the lines going through the same two points, 
but inclined at a certain angle to the vertical that is a little bit less than 7r/4. The 
results of EH and [2] are illustrated in [5j; see pp. 125, 126 and 137, respectively. 

In [3], another geometric condition alternative to (ii) was found which together 
with (i) guarantees the absence of trapped modes at all frequencies for fixed bodies. 
This condition does not impose any restriction on the angle between the surface¬ 
piercing contour and the free surface (arbitrarily small angles are admissible), but 
this is achieved at the expense that the wetted contour is subject to a certain point- 
wise restriction (it must be transversal to curves (1201) in a certain definite fashion). 

On the other hand, condition (i) is essential for the absence of trapped modes. 
This became clear when 8] constructed an example of such a mode for which pur¬ 
pose she applied the so-called semi-inverse method (see, for example, [7j for its brief 
description). Her example involves two fixed surface-piercing cylinders each of which 
satisfies the modified condition (ii) of EH, but they are separated by a nonzero spac¬ 
ing. Another example of a mode trapped by two fixed surface-piercing cylinders was 
found by ED]. Subsequently, [4] proved that the latter cylinders can be considered as 
two immersed parts of a single body which freely floats in trapped waves, but remains 
motionless. 

During the past decade, the problem of the coupled time-harmonic motion of water 
and a freely floating rigid body has attracted much attention. Along with the just 
mentioned paper [4], rigorous results were obtained in 17], where a brief review of 
related papers is given. However, the substantial part of work concerns the study of 
trapped modes and the corresponding trapping bodies and only the paper [6] has been 
focussed on conditions eliminating trapped modes in the case when a surface-piercing 
or totally submerged body is present (for a surface-piercing body the original proof 
of [2] was essentially simplified). In the present paper, our aim is to fill in this gap at 
least partially. 

In the present note, we find conditions on the frequency so that they guarantee 
that no modes (or some specific modes) are trapped by a freely floating body provided 
its geometry satisfies the assumptions used in EH and {3] for establishing the absence 
of modes trapped by the same body being fixed. 
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Figure 1: A definition sketch of the submerged cross-section of a cylinder. 

2 Statement of the problem 

Let the Cartesian coordinate system (x, y) in a plane orthogonal to the generators of a 
freely floating infinitely long cylinder be chosen so that the y-axis is directed upwards, 
whereas the mean free surface of water intersect this plane along the x-axis, and so 
the cross-section W of the water domain is a subset of R^ = {a; £ R, y < 0}. Let 
B denote the bounded two-dimensional domain whose closure is the cross-section a 
floating cylinder in its equilibrium position. Let both the immersed part B = BC\ R 2 . 
and the above-water part B \ Ri be nonempty domains and D = B C I SR 2 . be 
a nonempty interval of the x-axis, say {x £ (—a, a), y = 0} (see figured]). We 
suppose that W is either R 2 , \ B when water has infinite depth (see figure [T]) or 
{x £ R, —h < y < 0} \ B, where h > bo = swp^ x y ^ B |y|, when water has constant 
finite depth. We suppose that W is a Lipschitz domain, and so the unit normal 
n pointing to the exterior of W is defined almost everywhere on dW. Finally, by 
S = dB fl R'L and F = <9R?_ \D we denote the wetted contour and the free surface 
at rest, respectively; if water has finite depth, then H = {x £ R, y = —h} is the 
bottom’s cross-section. 

For describing the small-amplitude coupled motion of the system it is standard 
to apply the linear setting in which case the following first-order unknowns are used. 
The velocity potential <f >(x,y;t) and the vector-column q(t) describing the motion of 
body whose three components are as follows: 

• q\ and g 2 are the displacements of the centre of mass in the horizontal and vertical 
directions, respectively, from its rest position (x®,?/ 0 )); 

• <j 3 is the angle of rotation about the axis that goes through the centre of mass 
orthogonally to the (x,y)-plane (the angle is measured from the x- to y-axis). 

We omit relations governing the time-dependent behaviour (see details in 0). and 
turn directly to the time-harmonic oscillations of the system for which purpose we 
use the ansatz 

(<f>{x,y,t),q(t)) = Re{e~ lult (<p{x,y), i*)}, (1) 

where ui > 0 is the radian frequency, ip £ Hl oc (W ) is a complex-valued function and 
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X G C 3 . To be specific, we first assume that W is infinitely deep in which case the 
problem for (ip, %) is as follows: 


V 2 <p = 0 in W, 

(2) 

d y p — vip = 0 onF, where v = co 2 /g, 

(3) 

3 

d n p = to N J x N oXi) on 5 > 

1 

(4) 

V p — > 0 as y —> — 00 , 

(5) 

1 |«9i 3.1 <yC> — ivp\ 2 ds = o(l) as b —> 00 , 

VKn{ |x|=6} 

(6) 

to 2 Ex = — to / pNds + gKx- 
Js 

(7) 


Here V = (d x ,d y ) is the spatial gradient, g > 0 is the acceleration due to gravity 
that acts in the direction opposite to the y-axis; N = (Ni, IV 2 , N 3 ) T (the operation T 
transforms a vector-row into a vector-column and vice versa), where (JVi, A ^) 7 = n. 
N 3 = (a; — a 'S°\y — y^) x n and x stands for the vector product. In the equations 
of body’s motion 0 , the 3x3 matrices are as follows: 

(I M 0 0 \ /0 0 0 \ 

E = 0 I M 0 and K = 0 I D 1° . ( 8 ) 

VO 0 I») \o If I? x +I*) 

The positive elements of the mass/inertia matrix E are 


I M =Po 1 


J ^ p(x, y) dxdy and I™ = p 0 1 p(x, y) (x - a: (0) ) +(y - y (0) ) dxdy 


where p(x, y) > 0 is the density distribution within the body and po > 0 is the 
constant density of water. In the right-hand side of relation 0 , we have forces and 
their moments. In particular, the first term is due to the hydrodynamic pressure, 
whereas the second one is related to the buoyancy (see, for example, m ; the non¬ 
zero elements of the matrix K are 


I D = f dx > 0, I x = f (x — £^)d:E, 

Jd J D 

Ixx= [ ( x ~ a; ( 0 ) ) 2 da; > 0 , ly = [ (y - y {0) )dxdy. 

Jd Js 

Note that the matrix K is symmetric. 

In relations ©: © and 0, uj is a spectral parameter which is sought together 
with the eigenvector (ip, x)- Since W is a Lipschitz domain and p G Hl oc (W), relations 
© © are, as usual, understood in the sense of the following integral identity: 

/ S/pVip dxdy = v / ipipd x + uj , tpN T xds, (9) 

iff ■/ F Js 

which must hold for an arbitrary smooth having a compact support in W. Finally, 
relations © and © specify the behaviour of <p at infinity. The first of these means 
that the velocity field decays with depth, whereas the second one yields that the 
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potential given by formula o describes outgoing waves. This radiation condition is 
the same as in the water-wave problem for a fixed obstacle (see, for example, 0). 

The relations listed above must be augmented by the following conditions con¬ 
cerning the equilibrium position: 

• The mass of the displaced liquid is equal to that of the body: I M = f B dxdy 
(Archimedes’ law); • The centre of buoyancy lies on the same vertical line as the 
centre of mass: f B (x — x^) dxdy = 0; • The matrix K is positive semi-definite; 
moreover, the 2x2 matrix K' that stands in the lower right corner of K is positive 
definite (see HD- 

The last of these requirements yields the stability of the body’s equilibrium position, 
which follows from the results formulated, for example, by [Tj, §2.4. The stability 
is understood in the classical sense that an instantaneous, infinitesimal disturbance 
causes the position changes which remain infinitesimal, except for purely horizontal 
drift, for all subsequent times. 

In conclusion of this section, we note that relations (0 and 0 must be amended 
in the case when W has finite depth. Namely, the no flow condition 

d y v = 0 on H (10) 

replaces 0, whereas v must be changed to ko in 0, where ko is the unique positive 
root of kg tanh(fco/i) = v. 


3 Equipartition of energy, trapped modes and con¬ 
ditions guaranteeing their absence 

3.1 Equipartition of energy 

It is known (see, for example, m § 2 . 2 . 1 ]), that a potential, satisfying relations 0, 
©, © and 0, has the asymptotic representation at infinity of the same type as 
Green’s function. Namely, if W has infinite depth, then 

<p(x,y) = A±(y) e lu \ x \+r±(x,y), where |r±| 2 , |Vr±| = 0([x 2 +?/ 2 ] _1 ) as x 2 +y 2 —»• oo, 

( 11 ) 

and the following equality holds 

vf {\ A +(y)\ 2 + \A_{y)\ 2 ) dy = -Im [ Tpd n <pds. (12) 

J- oo Js 

Assuming that (<p, x) is a solution of problem © 0, we rearrange the last for¬ 
mula using the coupling conditions 0 and 0- First, transposing the complex con¬ 
jugate of equation 0, we get 

w 2 (ExV = -w [ TpN J ds + g (Kxf ■ 

Js 

This relation and condition 0 yield that the inner product of both sides with x can 
be written in the form: 

u 2 X J Ex - 9 X T Kx = - [ pd n ipds. (13) 

Js 
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Second, substituting this equality into (fl2l) . we obtain 

vj (\A + (y)\ 2 + \A_{y)\ 2 )dy = lm{uj 2 x T Ex-gX T Kxy (14) 

In the same way as in [2], this yields the following assertion about the kinetic and 
potential energy of the water motion. 

Proposition 1 . Let (<p,x) be a solution of problem ©-©, then 


/ |V(^| 2 dxdy < oo and v / \ip\ 2 dx < oo. 

Jw Jf 

Moreover, the following equality holds: 


(15) 


f \Vp\ 2 dxdy + uj 2 x T Ex = v f \p\ 2 dx + g x T Kx- (16) 

Jw J F 

Here the kinetic energy of the water/body system stands in the left-hand side, 
whereas we have the potential energy of this coupled motion in the right-hand side. 
Thus the last formula generalises the energy equipartition equality valid when a fixed 
body is immersed into water. Indeed, X = 0 for such a body, and (flbl) turns into the 
well-known equality (see, for example, formula (4.99) in [5]). 

Proposition 1 shows that if (<p, x) is a solution of problem © © with complex¬ 
valued components, then its real and imaginary parts separately satisfy this problem. 
This allows us to consider (<p, x) as an element of the real product space H 1 (W) x 1R 3 
in what follows (the sum of two quantities (11511 defines an equivalent norm in H 1 (W)). 


Definition 1. Let the subsidiary conditions concerning the equilibrium position 
(see § 2) hold for the freely floating body B. A non-trivial real solution (p, x) £ 
H l (W) x R 3 of problem © and © is called a mode trapped by this body, whereas 
the corresponding value of u> is referred to as a trapping frequency. 

In order to determine when (p, x) £ H l (W) x R 3 is not trapped by B we write 
m as follows: 


X T (w 2 Ex - g K)x = v [ \p\ 2 dx — [ |V^| 2 dxdy. (17) 

J F JW 

It is clear that the left-hand side is non-negative provided ui 2 is sufficiently large, and 
so we arrive at the following. 

Proposition 2. Let E and K be given by © and let w 2 be greater than or equal to 
the largest A satisfying det(A E — gK) = 0. If the domain W is such that the inequality 


v f l^l 2 da: < f |Vv?| 2 daxh/ (18) 

J f Jw 

holds for every non-trivial p £ H l {W), then uj is not a trapping frequency. 

Note that if W has finite depth, then v must be changed to fco in relation CD, 
where the behaviour of the remainder must be also replaced by the following one: 

|r±|, |Vr±| = 0(|a;|~ 1 ) as |a;| —» oo. (19) 

In relations (fl2l) and (fLll) v must be also changed to kg. On the other hand, formula 
CD remains be valid in the same form as above, and so proposition 2 is true in this 
case as well. 
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3.2 Examples of water domains for which inequality (fT8|) holds 

We begin with the case when W has infinite depth. By id and i-d we denote the 
rays emanating at the angle 7 r /4 to the vertical from the points (d, 0) and (— d, 0), 
respectively, and going to the right and left, respectively. 

Let the whole rays id and i-d belong to W for all d > a. Thus, B is confined 
within the angular domain between the lines inclined at 7 r /4 to the vertical and going 
through the points (a,0) and (—a, 0) to the right and left, respectively. Under this 
assumption, m proved (see also [5], §§3.2.2.1 and 3.2.2.2) that the inequality 

v / \p\ 2 dx < / |V(^| 2 dxd?/ 

Jf JWc 

holds provided ip satisfies conditions m and relations m and ©■ Here W c is the 
subset of W covered with rays {id : (d, 0) £f}U {i-d '■ (— d, 0) € F}. According 
to the last inequality, if ip is non-trivial, then (1181) holds. Therefore, proposition 2 is 
applicable, thus giving a criterion which values of w are non-trapping frequencies for 
the freely floating B whose immersed part B is confined as described above. 

In order to obtain inequality (fT51) in the case when W has finite depth, id and i-d 
must be replaced by similar segments connecting F and H and inclined at a certain 
angle to the vertical that is a little bit less than 7 r/ 4 . Numerical computations of [TTj 
show that the same result as for deep water is true when B is confined between the 
segments inclined at 441 . 

4 Another criterion eliminating some particular 
trapped modes 

In this section, we turn to the case when B does not satisfy the conditions of §3.2. 
To be specific, we suppose that W is bounded from below by the rigid bottom H. 
Moreover, we assume that B is symmetric about the y-axis (see figure [TJ) ; this implies 
that N\ = n x (N 2 = n y ) attains the opposite (the same, respectively) values at every 
pair of points on B which are symmetric about the y-axis. Let also p(x, y) be an even 
function of x, and so x ^ = 0 (the centre of mass lies on the y-axis); this implies that 
IV 3 = xn y — n x (y — y^) has the same behaviour as Ni. 

The last restriction on B or, more precisely, on B is expressed in terms of the 
curves 

x 2 + (y — acot tr ) 2 = a 2 (cot 2 a + 1 ), ±x > 0 , y < 0 , ( 20 ) 

parametrised by a £ (— 7 r, 0). On curves of these two families we define directions 
as shown in figure [T] It is clear that all curves (EU1) . that intersect H transversally, 
enter into W. Let this property also hold on S'; that is, all transversal intersections of 
curves (12U1) with S are points of entry into W (see figure [TJ. In what follows, a body 
satisfying the listed conditions is referred to as belonging to the class B provided the 
conditions considered in § 3.2 are not fulfilled for it. 

The following assertion generalises the criterion of [3] guaranteeing the absence of 
trapped modes for fixed surface-piercing bodies immersed in deep water and satisfying 
the above transversality condition with the family of curves (I2U1) . As in proposition 2 
the values of w that are not trapping frequencies must be sufficiently large, but what 
is new that some restrictions must be also imposed on the type of mode. 
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Proposition 3. Let W have finite depth and let B be a freely floating body belonging 
to the class B. If to 2 is strictly greater than the largest A such that det(A_E — gK) = 0 
with E and K given by © , then oj is not a trapping frequency for modes of the form: 

(a) <p is an even function of x and x — (c^i, 0 , c^ 3 ) T ; 

(b) ip is an odd function of x and x = ( 0 ,^ 2 , 0 ) T . 

Proof. Let us write relations © a and m using the bipolar coordinates (u,v). 
The corresponding conformal mapping is usually defined as follows (see, for example, 
@,§ 10 . 1 ): 


x = a sinh it/ (cosh it — cosu), y = asrnv / (cosYlu — cosv). ( 21 ) 

Therefore, m maps the strip {— oo < u < +oo, — n < v < 0} onto 1R?_ so that for 
every a € (— n, 0) the image of the left (right) half-line {±u > 0, v = a} is the circular 
arc (l20l) that lies in the left (right) half-plane (see figure©. Moreover, 

{—oo < u < + 00 , v = — 7r} and {±w > 0, v = 0} 

are mapped onto {|x| < a, y = 0} and {±x > a, y = 0} respectively. Finally, we have 
that 

|z , (£)| = a/(coshw — cosi>), where z = x + iy and £ = u + iv. 

The inverse mapping £(z) has the following properties: the points a and —a on the 
x-axis go to infinity on the ((-plane, whereas z = oo goes to £ = 0; thus F is mapped 
onto the whole it-axis. 

Denoting by W the image of W, we see that apart from the it-axis the boundary 
dW includes the images of S and H , say S and H respectively. According to properties 
of (Ell) , if B belongs to the class B, then S is symmetric about the x-axis, lies within 
the strip {—oo < u < +oo, —n < v < —a} and asymptotes the line v = —a as 
u — > ±oo; here a G ( 0 , 7 r) is the angle between S and F at (±a,0). Moreover, the 
right half of S is the graph of a decreasing function of u £ (0,+oo); its maximum 
value Vb is the root of cosv — (a/bo) sin v = 1- Finally, B is a closed curve with the 
following properties. It is symmetric about the x-axis, is tangent to the u-axis at the 
origin and is the graph of a concave function of v G (vh, 0 ); here Vh G (— 7 r, 0 ) is the 
root of cosx — ( a/d ) sinn = 1. It is clear that — 7 r /4 < —a < Vb < Vh < 0. 

Let f(u,v) = <p(x(u,v),y(u,v)), then relations © © yield that 

waNlx 

V (f = 0 in W, (coshu — l)<f> v = va <f> when v = 0, ■ rv = - - - -- on S. 

cosh u — cos v 

( 22 ) 

Here is the unit normal to S U H exterior with respect to W and = N z (q- 
Moreover, condition m implies that 


V 0 • n^- = 0 on H, 

whereas condition © takes the form 

(j)N(; dsf 


ui Ex = — w 


L 


s cosh u — cos V 


gK X - 


(23) 


(24) 
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Furthermore, conditions (lT5l) give that 



| V</)| 2 dudu < oo 


and 



cosh u— 1 


dw < oo, 


(25) 


whereas equality m turns into the following one: 

/ |V</>| 2 dwdn — va [ — } — d u = ~x T (<+ 2 E — g K)x- 

Jw J- oo cosh u — 1 


(26) 


Further considerations are based on the following identity (see 15], Subsection 

2 . 2 . 2 ): 

(2 u<t> u + <)>)V 2 0 = V • (2 u<t> u + d>)V</> - 2 eft - (u|V0| 2 ) u . (27) 

Here the left-hand side vanishes due to the Laplace equation for </>. Let us integrate 
this identity over W' = W n {|u| < b} and b is sufficiently large (in particular, 
H C (M < b}). Using the divergence theorem, we get 


2 f (j>^dudv+ ( u • n |V</>| 2 dS 

Jw 1 Js'UH 


r+b 


[2 u(j> u (u, 0) ± (j>(u, 0)] ip v (u, 0) du 


l-b 


+ / (2 u(j) u + (j>) V(/> • dS + ± / (2 u<j> u + 4>) <f> u du, (28) 

Js 1 ± Jc± 


where S' = S D (|u| < 5}, u = (zt, 0), denotes the summation of two terms 
corresponding to the upper and lower signs, respectively, and C± = W D (zt = ±6}. 
All integrals on the right arise from the first term on the right in G71) and one more 
integral of the same type vanishes in view of the boundary condition (1231) on "H. 

Let us consider each integral standing on the right in (E51) . Using the free-surface 
boundary condition, we get that the first term is equal to 


p-\-b 

va / [2u(/) u (u, 0) + 4>{u, 0)] 

J-b 


&{u, 0 ) 

cosh u—1 


du 


r+b 


. usinhud> 2 (u, 0) , 

= m h (c.sh„-l)^ d “ + TO 


u </> 2 (u, 0) 


cosh u—1 


L L= — b 


where the last expression is obtained by integration by parts. It follows from (1151) that 
<p(x, y) tends to constants as ( x , y) —» (±a, 0), and so <j>(u, v) has the same property as 
u -+ ±oo. Therefore, the integrated term in the last equality tends to zero as b —> oo, 
whereas the integral on the right converges in view of (1251) . 

The second integral on the right in (11251) is equal to 

f (2u0„ + <t>)Njx . „ 

toa / -^- 2 —do. 

J S' cosh u — cos V 

Since S belongs to the class 6, we have that <j> and ip are simultaneously even and 
odd functions of x and u respectively. Therefore, either of the assumptions (a) and 
(b) implies that this integral vanishes because the integrand attains opposite values 
at points of S' that are symmetric about the u-axis. 

Finally, (l25l) implies that there exists a sequence tending to the positive 

infinity and such that the last sum in (l28l) tends to zero as bk —> oo. Passing to the 
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limit as k —> oo, we see that the transformed equation ( 1281 ) with b = bk gives the 
following integral identity: 

2 [ (fadudv + [ u • n |V0| 2 dS” — va [ - —0) du = 0 

Jw Jsun J- oo (coshu-1) 2 

provided either of the assumptions (a) and (b) holds. 

Subtracting this from ( 1261 ) multiplied by two, we get 


f (ftldudv— f u • n |V</>| 2 d5 
Jw Jsun 


/w 

c+oo 


+ua 


:sinhu — 2(coshu — 1) 
(coshu — l) 2 

= -2 X T (c j 2 E- 9 K) X . 


4> 2 {u , 0) du 


(29) 


If u} 2 is strictly greater than the largest A such that det(A.E — gK) = 0, then cannot 
hold unless lo is not a trapping frequency for modes of the form (a) and (b). Indeed, 
the right-hand side is negative for such a value of u> and a non-trivial Xi whereas the 
left-hand side is non-negative because S belongs to the class B and the fraction in the 
last integral is non-negative. The obtained contradiction proves the proposition. □ 


5 Conjecture 

Given the proof of a theorem guaranteeing the uniqueness of a solution to the lin¬ 
earised problem about time-harmonic water waves in the presence of a fixed obstacle, 
then this proof admits amendments transforming it into the proof of an analogous 
theorem for the same obstacle floating freely with additional restrictions on the non¬ 
trapping frequencies (they must be sufficiently large) and, in some cases, on body’s 
geometry and on the type of non-trapping modes. 
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